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We extend our approach for the exact solution of the Riemann problem in relativistic hydrody- 
namics to the case in which the fluid velocity has components tangential to the initial discontinu- 
ity. As in one-dimensional flows, we here show that the wave-pattern produced in a multidimen- 
sional relativistic Riemann problem can be predicted entirely by examining the initial conditions. 
Our method is logically very simple and allows for a numerical implementation of an exact Rie- 
mann solver which is both straightforward and computationally efficient. The simplicity of the 
approach is also important for revealing special relativistic effects responsible for a smooth tran- 
sition from one wave-pattern to another when the tangential velocities in the initial states are 
suitably varied. While the content of this paper is focussed on a flat spacetime, the local Lorentz 
invariance allows its use also in fully general relativistic calculations. 



1. Introduction 

Since the solution of the Riemann problem has been introduced in numerical hydrodynam- 
ics (Godunov 1959), high resolution shock capturing methods have become a common tool for 
handling nonlinear hydrodynamical waves in a variety of physical situations and the subject of 
detailed mathematical analyses both in Newtonian (Leveque 1992, Toro 1997) and relativistic 
regimes (Smoller & Temple 1993). Nonlinear waves of this type are very common in astro- 
physical scenarios (such as gamma-ray bursts, accretion onto compact objects, relativistic jets, 
supernova explosions) in which the motion of the fluid is characterised by relativistic speeds and 
by the appearance of strong discontinuities (see Ibanez & Marti, 1999 and Marti & Miiller, 1999 
for a review). In a recent paper, Rezzolla & Zanotti (2001, hereafter paper I) have proposed a new 
procedure for the exact solution of the Riemann problem which uses the relativistically invariant 
relative velocity between the unperturbed "left" and "right" states of the fluid in order to extract 
important information contained in the initial data and that more traditional approaches were not 
able to put into evidence. Most notably, it was shown that, given a Riemann problem with as- 
signed initial conditions, it is possible to determine in advance both the wave-pattern that will be 
produced after the removal of the initial planar discontinuity and the bracketing interval of the un- 
known pressure in the region that forms behind the wave fronts. Besides clarifying some aspects 
of the Riemann problem, the approach proposed by Rezzolla & Zanotti, and that was foreseen by 
Landau & Lifshitz (1959) in Newtonian hydrodynamics, was shown to be computationally more 
efficient as compared to the more traditional ones. 

In a recent work, Pons et al (2000) have extended the analytic solution of the one-dimensional 
Riemann problem (Marti & Miiller 1994) to the case in which nonzero velocities tangential to 
the initial planar discontinuity are present. An important result obtained by Pons et al (2000) 
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was to show that the introduction of tangential velocities and the appearance of global Lorentz 
factors linking quantities on either side of the discontinuity can affect the solution of the Riemann 
problem considerably. In the present paper we show that the approach by Rezzolla & Zanotti can 
be successfully extended to this more general case and with the same advantages that were found 
in the case of zero tangential velocities. Within this new procedure, the introduction of tangential 
velocities is imprinted in the expression for the jump of the velocity normal to the discontinuity 
surface, and this has allowed us to reveal interesting special relativistic effects. In relativistic 
hydrodynamics, in fact, the wave-pattern produced by the decay of the initial discontinuity can 
be changed by simply varying the tangential velocities in the initial states, while keeping the 
rest of the thermodynamic quantities of the Riemann problem unmodified. This effect has no 
analogue in Newtonian hydrodynamics (Rezzolla & Zanotti 2002). 

The plan of the paper is as follows: after a review of the method in Section 2, we report in 
Section 3 the hydrodynamical equations relevant for the present discussion of nonzero tangential 
velocities. In Section 4 we show how to use the invariant relative velocity to extract from the 
initial data the information on the wave-pattern produced, while Section 5 is devoted to the 
presentation of the special relativistic effects. The conclusions of the paper are in Section 6, 
and two Appendices complete the discussion providing the mathematical details of the results 
obtained in the main text. 

We use in this paper a system of units in which c = 1, a space-like signature (— , +, +, +) and 
a Cartesian coordinate system (t, x,y, z). Greek indices are taken to run from to 3, Latin indices 
from 1 to 3 and we adopt the standard convention for the summation over repeated indices. 



2. A brief review of the method 

In a flat spacetime consider a perfect fluid described by the stress-energy tensor: 

T» v ee (e + p)u^u u + prf v = phu^u u + prf v , 



(2.1) 



where rf iV = diag(— 1, 1, 1, 1) and e, p, p, and h are the proper energy density, the isotropic 
pressure, the proper rest mass density, and the specific enthalpy, respectively. 

Further assume the fluid to consist of an initial "left" state (indicated with an index 1) and 
an initial "right" state (indicated with an index 2), each having prescribed and different values 
of uniform pressure, rest-mass density and velocity. The two discontinuous states are initially 
separated by a planar surface S placed at a constant value of the x coordinate so that the unit 
space-like 4-vector ?t.q normal to this surface at t = has components = (0, 1, 0, 0)[{| Because 
we are considering a multidimensional flow, the fluid 4-velocity on either side of the initial 
discontinuity is allowed to have components in spatial directions orthogonal to n,Q, i.e. 



u" ee W(l,v x ,v y ,v z ) , 
J 1 is the square of the Lorentz factor and v 2 



(2.2) 



where W^ 2 = (1— u 2 ) -1 is the square of the Lorentz factor and v 2 = v l Vi = (v x ) 2 + (v v ) 2 + (v z ) 2 
is the norm of the 3-velocity. The hydrodynamical properties of the initial left and right states are 
described by the "state-vectors" 



Ui, a = 



I P \ 



where we have indicated with v* = [(v v ) 2 + (v 2 ) 2 ] 1 / 2 the tangential component of the three 

f More precisely, the unit normal to the hypersurface So is defined as the one-form mapping each vector 
tangent to the surface into zero. 
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velocity, satisfying the obvious relativistic constraint that (i>*) 2 + (v x ) 2 < 1 (Note that in paper 
I the 4-velocity vector was chosen so that v v = v z = = u*). Hereafter, we will refer to v x as 
the normal velocity. 

The fluid states Ui and U2 represent the initial conditions of a multidimensional '"Riemann 
problem" whose solution consists of determining the flow that develops when the system is 
allowed to relax. In general, the temporal evolution can be indicated as (Marti & Miiller 1994) 

LW^L*CR*W^R , (2.3) 

where W denotes a nonlinear wave (either a shock, S, or a rarefaction wave, TV), propagating 
towards the left (<— ) or the right (— >) with respect to the initial discontinuity. Moreover, L„ and 
R* are the new hydrodynamic states that form behind the two nonlinear waves propagating in 
opposite directions. A contact discontinuity, C, separates the region L* and and it is charac- 
terised by the fact that both the pressure and the normal velocity are continuous across it, while 
both the rest-mass density and the tangential velocities can be discontinuous. 

The approach introduced in paper I focuses on (vf 2 )o, the relativistic invariant expression for 
the relative velocity between the two unperturbed initial states. By construction, this quantity 
measures the relativistic jump of the fluid velocity normal to the discontinuity surface. The so- 
lution of the relativistic Riemann problem is then found after the pressure in the region between 
the two nonlinear waves, p*, is calculated as the root of the nonlinear equation 

»i 2 (P.) - « 2 )o = , (2.4) 

where vf 2 (p*) has a functional form that is different for each of the three possible wave -patterns 
that might result from the decay of the initial discontinuity. The key aspect of the new approach is 
that the wave-pattern produced by the decay of the initial discontinuity can be entirely predicted 
in terms of the initial data U1.2. This represents an important advantage since it allows to deduce 
in advance which set of equations to use for the solution of the exact Riemann problem and the 



interval bracketing the root of (2.4) 



The validity of the approach discussed in paper I is in the mathematical proof that the func- 
tion vf 2 = v f 2 (p*) i s monotonically increasing with and it is composed of three branches 
corresponding to the three possible wave-patterns. Furthermore, the three different branches al- 
ways join smoothly through specific values of vf 2 (p*), denoted respectively as (vf 2 ) 2S , (vf 2 ) SR , 
(v x 2 ) 2 J^(cf. Fig. 1 of paper I). In view of these properties, it is possible to compare (vf 2 )o with 
the relevant limiting values (vf 2 ) 2S , (vf 2 ) SH1 (i 'f 2 ) 2R constructed from the initial conditions and 



determine, prior to the solution of equation ( 2.4 ), both the wave-pattern that will be produced and 



the functional form of vf 2 (p» ) to be used. The logical scheme to be followed in this procedure is 
synthesised as follows 



(vf 2 ) > (vf 2 ) 2S ■ LS^LXR*S^R "12 (p*) = « 2 ) 2S 

(vf 2 ) S R < K 2 )o < (612)25 : LK^L*CR*S^R, vf 2 (P*) = 0i 2 ) SH 

< « 2 )o < («ia) M : LU^L*CR*H^R, »f 2 (p.) = (vf 2 ) 2R 

(^12)0 < (vi 2 ) 2R '■ LIZ^L^CR^IZ^R with vacuum, — 



(2.5) 

After the initial relative normal-velocity {vf 2 )o has been compared with the limiting values for 
the different branches and the correct wave-pattern [i.e. the functional form v^Cp*)] has been 
found, the numerical solution of the exact Riemann problem can be performed. Furthermore, 
the knowledge of the wave-pattern relevant for the solution of the problem provides the values 

f Here 2S, SR and 2R indicate a two shock, a shock-rarefaction and a two rarefaction wave-pattern, 
respectively. 
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p max and p m in bracketing the unknown pressure p*. In the case of a one-dimensional flow, the 
bracketing interval is sinthesized in the scheme below 





LS^- LifCR* 5_> R 


LlZ^L^CRxS^R 


LTZ^LXR*n^R 




max(pi,p 2 ) 


min(pi,p 2 ) 





Pmax 


oo 


max(pi,p 2 ) 


min(pi,p 2 ) 



In practice, this novel approach has two major advantages. Firstly, it results into a straightforward 
numerical implementation and in a more efficient numerical algorithm with a reduction of the 
computational costs. Secondly, it simplifies considerably the logic formulation of the relativistic 
Riemann problem. As we shall discuss in Section ||, this feature has been essential in revealing 
new special relativistic effects. 

The extension of the approach presen ted i n paper I to the case when tang enti al velocities are 
present is straightforward and equation (2.4) as well as the logical scheme (2.5) apply unmodi- 
fiedf . The only changes introduced by the presence of tangential velocities are restricted to the 
expressions for the limiting values of the relative velocity (vf 2 ) 2S , {v'f 2 ) SR , {v^ 2 ) 2R . The de- 
tails of these changes will be presented in the following Sections which have been written for a 
generic equation of state (EOS) and use an ideal fluid EOS as a test case. 



3. Hydrodynamical Relations across the waves 

As discussed in the previous Section, the expression for the relative normal-velocity between 
the two initial states of the Riemann problem represents the building block in our approach and, 
to simplify our notation, hereafter we will refer to the different flow regions using the following 
mapping 

LW^L£R*W^R 1 3 C 3' 2 . (3.1) 

While the values of vf 2 are relativistic invariants, there exists a reference frame which is better 
suited to evaluate this quantity. In the reference frame of the contact discontinuity, in fact, the 
normal velocities behind the nonlinear waves are, by definition, zero (i.e. v'f c — = uf, c ) and 
the relative velocities across the nonlinear waves measured in this reference frame will be 

Because of their invariance, the normal velocity jumps across the nonlinear waves measured in 
the Eulerian frame can be expressed as 

f Note that in the case of two shock waves propagating in opposite directions the maximum 
value for the pressure needs not be infinite but can be determined from the solution of the equation 
vfzip*) ~ ("?2)max = 0, where (vf 2 )max is the maximum allowed value of the normal relative veloc- 
ity that is compatible with the assigned initial tangential velocities. 
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fe) = ^J = fe),c=^, C . (3.5) 
i v 2 u 3 , 

As a result, the relative normal-velocity between the two initial states can be written as 



,J i,C - V 2,C 

i - Kc)Kc) 



= = r a s ■ 0.6) 



In what follows we will briefly discuss how to calculate the normal-velocity jump across a 
shock wave and a rarefaction wave, respectively. The expressions derived in this way will then 
be used to calculate vf c and v% c necessary to build vf 2 = vf 2 (p*) [cf. eq. (j3~6|)]. 

3.1. Jumps Across a Shock Wave 

Following Pons et al (2000), we use the Rankine-Hugoniot conditions in the fixed Eulerian ref- 
erence frame to express the normal velocity of the fluid on the back of the shock front in terms 
of the pressure^ 

x = h a W a v* + W s (p b - Pa )/J 

b KW a + (p b -p a )[W s vZ/J+l/(p a W a )} ' 
where J is mass flux and the subscripts b and a denote a quantity evaluated on the back or ahead 
of the wave front, respectively. In equation (3.7), W s = (1 — V^ 2 ) -1 / 2 is the Lorentz factor of 
the shock velocity V s , with the latter being 



and with the ± signs referring to a shock wave propagating to the right or to the left, respectively. 



Besides giving the jump in the normal velocity across a Shockwave, expression (3.7) states that 
the two regions of the flow across the shock wave are effectively coupled through a Lorentz factor 
which, we recall, is built also in terms of the tangential velocities. This is a purely relativistic 
feature and an important difference from Newtonian hydrodynamics, in which the solution of 
the Riemann problem does not depend on the tangential component of the flow. Many of the 
consequences introduced by this coupling will be further discussed in Section || but one of them 
can already be appreciated when looking at the jumps across the shock wave of the tangential 
velocities that can be written as 

[[hWv y ]} = = [[hWv z ]} , (3.9) 
and where we have adopted the standard notation in which [[F]] = F a — Ft, (Anile 1989). 



Conditions (3.9 (basically state that the ratio v v jv z remains unchanged through rarefaction waves 
and therefore that the tangential velocity 3-vector does not rotate but can change its norm. This 
property, which applies also across rarefaction waves, represents a major difference from the 
behaviour of the tangential velocities across Newtonian nonlinear waves. In this latter case, in 
fact, the tangential velocity vector does not change its norm, nor rotates: [[v v ]] = = [[v z ]]. 
The square of the mass flux across the wave can be expressed as 

a M] 

J =-wkv OM} 

f Note that when tangential velocities are present, calculating jump conditions in the rest frame of the 
shock front is not particularly advantageous. In this case, in fact, the velocity jump across the shock cannot 
be expressed as an algebraic relation among the thermodynamical quantities across the shock [cf. eq. (3.1) 
of Rezzolla & Zanotti 2001]. Rather, the ratio of the velocities ahead and behind the shock front need 
to be found as root of a nonlinear equation. Analytic solutions to this equation can be found only in the 
weak-shock limit or for an ultrarelativistic equation of state (see Koenigl 1980 for a discussion). 
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M. 



(3.11) 



where the ratio h/p in the shocked region can be calculated through the Taub adiabat (Taub, 
1948) 

'ha hb 
.Pa Pb, 

In a general case, the mass flux can be obtained as a function of just one thermodynamical 
variable (p*) after using the EOS and the physical (h > 1) solution of the nonlinear equation 
(3.11). In the case of an ideal fluid EOS, 

p = (7 - l)pe = k{s)pi, (3.12) 

where 7 is the adiabatic index, and k(s) is the polytropic constant (dependent only on the specific 
entropy s), this can be done explicitly because ( 3 - 1 Oj ) and (3.11) take respectively the simple form 



j = — 



7-1 [[&(*- i)/p]] ' 



and 



(7- l)(Po - Pb) 



hi 



(7- !)(Pa -Pb) , . h a (p a -p b ) 
—rib + - 



(3.13) 



(3.14) 



IPb J " IPb Pa 

3.2. Jumps Across a Rarefaction Wave 

When considering a rarefaction wave it is convenient to introduce the self-similar variable 
£ = x/t in terms of which similarity solutions to the hydrodynamical equations can be found. 
An explicit expression for £ can be obtained by requiring that non-trivial similarity solutions for 
the rarefaction wave exist. This then yields (see Pons et al 2000 for details) 



- C*) ± C s ^(l-V 2 )[l~V 2 C 2 s -(v x ) 2 (l-C 2 )} 



1 



(3.15) 



where here too the ± signs correspond to rarefaction waves propagating to right or to the left, 
respectively. In the case of a perfect fluid, the isentropic character of the flow allows to calculate 
the velocity on the back of the wave as a solution of an ordinary differential equation 

dv x _ 1 (1 -& x ) 
dp 



(3.16) 



phW 2 (£-u x ) 

In principle, to calculate the normal fluid velocity at the tail of the rarefaction wave one should 
solve the ordinary differential equation (3.16), which might be very expensive numerically. To 
overcome this, it is convenient to make use of constraints such as those in expressions (3.9) 
(which remain valid also across a rarefaction wave) and express equation ( 3.16| ) in a different 
way. Defining A = h a W a v a — hbWbv\, the tangential velocity along a rarefaction wave can be 
expressed as 

'1 - (v x ) 2 ~ 



,.t\2 



A 4 



h 2 + A 2 



(3.17) 



This allows us to eliminate the depen dence on v t from equation ( 3.15 ). From the definition of 
the Lorentz factor and equation (3.17) it is straightforward to obtain 



W 2 = 



h 2 + A 2 



h 2 (l-(v x ) 2 ) ' 



and after some algebra one can arrive at 

(l-fo 35 ) 
(€ " v x ) 



= ± 



^h 2 +A 2 (l-c 2 ) 
hc s 



(3.18) 



(3.19) 
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Using this results, equation (3.16) can be written as follows: 



dv x y//l 2 +^ 2 (l- C g) dp 

l-(^) 2 {h 2 + A 2 ) pc s ' } 

Note that in this way we have isolated the thermodynamical quantities on the right hand side of 
( 3.20| ) and the kinematical ones on the left hand side, which can then be integrated analytically. 



For some particular cases (for example when the sound speed is constant), the right hand side 
too is integrable but for a generic EOS a numerical integration is necessary. The velocity at the 
tail of the rarefaction wave can then be obtained directly as 



where 



u£ = tanhS, (3.21) 

Here, h = h(p, s), p = p(p, s), and c s = c s (p, s), and the isentropic character of rarefaction 
waves allows to fix s = s a . Despite its complicated look, the integrand is a smooth, monotonic 
function of p, and a Gaussian quadrature with only 10-20 points has proved to be more accurate 
and efficient than a third order Runge-Kutta integrator requiring hundreds of function evaluations 



to solve ( |3.16| ). Testing and details on the numerical implementation of this procedure will be 
reported in a separate work. 

A couple remarks should now be made about the flow properties of a multidimensional and 
relativistic Riemann problem. The first remark is about the flow conditions at the contact dis- 
continuity where, as mentioned before, not only the rest-mass density but also the tangential 
velocities can be discontinuous. It is well known that under these conditions a Kelvin-Helmholtz 
(or shear-layer) instability could develop (Chandrasekhar, 1961), reach the nonlinear waves if 
these move subsonically and thus destroy the solution of the Riemann problem. While this could 
be a concern in principle, it is never a problem in practice. The reason for this is simple and is 
related to the fact that any implementation of an exact Riemann solver in a numerical code will 
provide the solution of the Riemann problem at a new time-level which is "Courant-Friedrichs- 
Lewy limited" from the previous one (Press et al 1986). Besides providing numerical stability, 
this time-step limitation also guarantees that the solution at the new time-level is everywhere 
causally connected with the one at the previous one. As a result, even if a Kelvin-Helmholtz in- 
stability would develop at the contact discontinuity propagating outward in the flow at the sound 
speed, this would not be able to influence the dynamics of the nonlinear waves over the time of 
interest in numerical calculations. 

The second remark is about the presence of tangential velocities across a shock-wave which 
can be in general discontinuous. While the stability of relativistic planar shocks in one-dimensional 
flows has been the subject of several investigations in the past (Anile & Russo 1987a, 1987b, 
Anile 1989), very little if anything at all is known about the stability properties of shock waves 
with discontinuous tangential velocities. It is reasonable to expect that a nonzero mass flux across 
these discontinuities would stabilise them against the Kelvin-Helmholtz instability, but further 
work needs to be done in this direction before a conclusion can be drawn. 



4. Limiting relative velocities 

As mentioned in the previous Sections, the basic operation in our approach consists of calcu- 
lating the relative normal-velocity across the two initial states and comparing it with the limiting 
relative velocities for each of the three possible wave-patterns. In practice, this amounts to calcu- 
lating equation (3.6) making use of expressions (3.4) and (3.5). In the following we will briefly 
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discuss the guidelines for the evaluation of the limiting relative velocities. In doing so we will 
adopt the convention of paper I and assume that p\ > p2, with the x— axis normal to the discon- 
tinuity surface being positively oriented from 1 to 2. 



4.1. lS^'iC'i' 2: Two Shock Waves 

We first consider a wave-pattern in which two shock waves propagate in opposite directions. In 
this case, the general expression for the relative normal-velocities between the two initial states 
(vf 2 ) 2S can be calculated from (3.6) wit h th e velocities behind the shock waves uf and Vy being 
determined through the jump condition (3.7). Because p\ is the smallest value that the pressure at 
the contact discontinuity p^ can take, the limiting value for the two shock waves branch (vf 2 ) 2S 
can be expressed as 



(S? 2 ) a 



lim [v\ 



(4.1) 



P3^Pl 



Evaluating the limit (4. 1 ) basically involves calculating the limits of v\ c and v 2 c for p^ tending 
to p\. Both these limits are straightforward to calculate and are 



lim r 

P3— >Pl 

lim v. 

P3^Pl 



l.C 



2.C 



= 0, 



to 



where Vy is simply the value of v% for p 3 



1 

Pi, i.e. 



V 2 V 3' 



vZ = lim v%, 

P3~*Pl 



(4.2) 
(4.3) 

(4.4) 



Using now the limits (4.2)-(4.3) and some lengthy but straightforward algebra, the explicit 
analytic expression for the limiting value of the two shock waves branch can be calculated as 



(3? 2 ) a 



lim v. 

P3^P1 



( P1 -P2){1-V* 2 V S ) 



2,C 



(V s - v$){h 2P2 (W 2 )*[l - {v x 2 Y)+ Pl - P2 } 



(4.5) 



Here V s is the velocity of the shock wave propagating towards the right in the limit of p$ — > 
pi and an explicit expression for it can be found in Appendix B in the case of an ideal fluid. 
Expression ( |4.5[ ) will be discussed further in Section || but it sufficient to point out here that the 
threshold value (vf 2 ) 2 s does not depend on the initial velocity in the state 1, v±. 



4.2. 1 1Z^- 3 C 3' 2: One Shock and one Rarefaction Wave 

We next consider the wave-pattern consisting of a rarefaction wave propagating towards the left 
and of a shock wave propagating towards the right. Also in this case , (vf 2 ) SR can be calculated 
from (3.6) with vf, bein g dete rmined through the jump condition (3.7) and from the numerical 
integration of equation ( |3.16| ) in the range [pi,P3]. Because p2 is now the lowest pressure in the 
unknown region behind the two nonlinear waves, the limiting value for the one shock and one 
rarefaction waves branch (vf 2 ) SR can be expressed as 



P3— >P2 

In the limit p% — > p 2 , the right-propagating shock is suppressed, uf, 



v 2 so that 



lim vie = , 



and 



P3^P2 



K 2 ) S h = lim vie . 

P3^P2 



(4.6) 



(4.7) 



(4.8) 
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Defining now 



Bi = I log ( \±4 ) , (4-9) 



2 ° V 1 



and using (3.21 ), it is readily obtained that 



/ m .lya + A 2 (l- c 2 ) \ 
(v x l2 ) SR = Km tanh(Bi - B) = tanh / V \\ s -±dp , (4.10) 

where the above integral can be evaluated numerically. A better look at the integral shows that 
only quantities in the left state are involved (through the constant A\ = h{W\v\) and that 
(vf 2 ) SR d° es n °t depend on the initial velocity in the state 2, v 2 . This property has an important 
consequence that will be discussed in Section || 

4.3. 1 7£<_ 3 C 3' 7£_> 2: Two Rarefaction Waves 
When the wave-pattern consists of two rarefaction waves propagating in opposite directions, 
(vf 2 ) 2R can De calculated from (3.6) with the velocities behind the waves being calculated using 
( fml ) and ( ^|| ). Since the lowest value of the pressure behind the tails of the rarefaction waves 



is zero, the limiting value for the two rarefaction waves branch (vf 2 ) 2R is given by 

(vf 2 ) 2R = lim (vf 2 ) 2R . (4.11) 

Proceeding as in previous subsection, we can now express {vf 2 ) 2R as 

r~x \ "i>c — v 2 ,c 



(«?,c)(«f,c) ' 



(4.12) 



where 



h i/ 0i f^iF*)' (4 - i3> 



y/fe 2 +^l(l 

{h 2 +Al)pc 



Ac = tanh v ^ 2 \_ _ - dp | , (4.14) 




and where A 2 = h 2 W 2 v\. While , the determination of {vf 2 ) 2R requires the numerical calculation 



of the integrals ( |4.13[ ) and ( |4.14| ), it has very little practical importance as it marks the transition 
to a wave-pattern with two rarefaction waves separated by vacuum; this is a very rare physical 
configuration which cannot be handled by a generic numerical code. 



Note that in computing ( XA2 ), both the left state quantities and the right ones are involved 
and, as a result, {vf 2 ) 2R will depend on both vi and v 2 . This property will be important in the 
subsequent discussion in Section ||. 

Fig. [j] shows the functional behaviour of vf 2 = vf 2 (p^) and how this behaviour is changed by 
the presence of nonzero tangential velocities. The initial conditions are those of a modified Sod's 
problem (Sod, 1978) in which Pl = 1.0, p x = 1.0, vf = 0.0, p 2 = 0.1, p 2 = 0.125, uf = 0.0, 
and 7 = 5/3. Each of the two curves shown is effectively the composition of three different 
ones (joined at the solid dots) corresponding to wave-patterns consisting of two shock waves 
(25*), one shock and one rarefaction wave {SR), and two rarefaction waves {2R). While the 
solid curve refers to initial conditions with zero tangential velocities, the dashed one is produced 
when nonzero tangential velocities, v\ — 0.7 = v 2 , are considered. Note that also in this latter 
case, the three branches are monotonically increasing with p% (a fundamental property whose 
mathematical proof can be found in Appendix A) but are all altered by the presence of nonzero 
tangential velocities. The consequences of this will be discussed in the next Section. 
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FIGURE 1. Relative normal-velocity between the two initial states as a function of the pressure at the 
contact discontinuity. Each curves is the continuous joining (marked by the solid dots) of three different 
curves corresponding respectively to two shock waves (25*), one shock and one rarefaction wave (SR), and 
two rarefaction waves (2R). The solid line refers to the case of zero tangential velocities, while the dashed 
line to the case in which v\ = 0.7 = v\. The initial state vectors are those of Sod's problem. 

5. Relativistic Effects 

The changes in the functional behaviour of vf 2 = vf 2 (P3) introduced by nonzero initial tan- 
gential velocities suggest that new qualitative differences could be found in a relativistic multi- 
dimensional Riemann problem. This was first discussed in a recent paper (Rezzolla & Zanotti 
2002) where the basic features of new relativistic effects were briefly pointed out. This section 
is dedicated to a more detailed discussion of how the changes in the functional behaviour of 
v i2 = v i2 (Ps) are responsible for relativistic effects in the dynamics of nonlinear waves. Before 
entering in the heart of the discussion, however, it is useful to remind that in Newtonian hydro- 
dynamics a multidimensional Riemann problem does not depend on the values of the tangential 
velocities in the two initial states. Rather, different wave-patterns can be produced only after a 
suitable change in either the normal velocity, the rest-mass density or the pressure. This is es- 
sentially due to the fact that tangential velocities are not changed across Newtonian nonlinear 
waves. In relativistic hydrodynamics, on the other hand, this is not the case and is at the origin of 
the effects discussed below. 

Let us restrict our attention to a situation in which the tangential velocity of only one of the 
two initial states is varied. This is simpler than the general case as it basically represents a one- 
dimensional cross-section of the three-dimensional parameter space, but it maintains all of the 
relevant properties. The two panels of Fig. || show the relative normal-velocity for the same initial 
conditions of Fig. [j] where either v\ or v\ is varied while all the other quantities of the initial state 
vectors are left unchanged. Different line types mark the different branches (joined at the filled 
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FIGURE 2. The same as in Fig. [IJbut here the different line types mark the different branches corresponding 
to two shock waves (dashed line), one shock and one rarefaction wave (dotted line), and two rarefaction 
waves (continuous line), respectively. The two panels show how the functional behaviour is modified when 
only one of the initial tangential velocities is varied (v\ for the left panel and v\ for the right one) while all 
the other components of the initial state vectors are left unchanged. 



dots) describing the relative velocity corresponding to two shock waves (25, dashed line), one 
shock and one rarefaction wave (SR, dotted line), and two rarefaction waves (2R, continuous 
line), respectively. Both panels of Fig. ^ indicate that when tangential velocities are present the 
relative normal-velocity is a function of p^ but also of v\ and i^IR 




0.2 0.4 0.6 0.8 



FIGURE 3. Transition from a SR wave-pattern to a 2S one. The left and right panels show the exact solution 
of the Riemann problem corresponding to models (a) and (e) in Table I, respectively. The initial state vectors 
are identical except for the values of v\. Solid, dotted and dashed profiles refer to pressure, rest-mass density 
and normal velocity, respectively. 

t For the of initial conditions chosen in Fig. ^| the position of (via) 2 n is very close to the limit — 1. 
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Consider, for instance, the case in which the normal velocities are chosen to be v'f = 0.5, wf = 
0.0, and that there are no tangential velocities. In this case, (vf 2 )o — 0.5 and the left panel of 
Fig. H shows that the solution to the Riemann problem falls in the SR branch, hence producing a 
wave-pattern consisting of a shock and a rarefaction wave moving in opposite directions. This is 
shown in more detail in the left panel of Fig. |^ where the different types of line show the solution 
of the Riemann problem at a time t > for the pressure (continuous line), the rest-mass density 
(dotted line) and the velocity (dashed line). 

However, if we now maintain the same initial conditions but allow for nonzero tangential 
velocities in state 2, the left panel of Fig. ^| also shows that the solution to the Riemann problem 
can fall in the 25 branch, hence producing a wave-pattern consisting of two shock waves moving 
in opposite directions. This is shown in the right panel of Fig. [| which illustrates the solution of 
the same Riemann problem but with initial tangential velocities v\ = 0.0 and v\ = 0.9. Note 
that except for the tangential velocities, the solutions in Figs. [|have the same initial state-vectors 
but different intermediate ones (i.e. p3, p^, p3',andwf). 

The Riemann problem shown in Fig. ^]is only one possible example but shows that a change in 
the tangential velocities can produce a smooth transition from one wave-pattern to another while 
maintaining the initial states unmodified. Furthermore, because the coupling among the different 
states is produced by the Lorentz factors this is not sensitive on the sign chosen for the tangential 
velocity. The transition from one wave-pattern to the other is better illustrated in Fig. [| where we 




FIGURE 4. Sequence of solutions for the pressure in Sod's problem. The initial tangential velocity I'Ljs 
gradually increased from to 0.9. The first and last solution of this sequence are also plotted in Fig. Q 

have collected in a three-dimensional plot a sequence of solutions for the pressure in which v\ is 
gradually increased from to 0.9. Note that when v\ = 0, the SR wave-pattern is well defined 
and the pressure at the contact discontinuity is intermediate between p\ and pi. Note also that as 
v\ is increased, the wave-pattern gradually changes, pa increases up until it becomes larger than 
Pi, signalling the transition to a 25 wave-pattern. 

Interestingly, the transition does not need to always produce a solution consisting of two shock 
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waves. Suppose, in fact, that the normal velocities are now chosen to be vf = 0.0, = 0-5- 
We can repeat the considerations made above and start by examining the wave-pattern produced 
when there are zero tangential velocities. In this new setup, (vf 2 )o — —0.5 and the right panel 




0.2 0.4 0.6 0.8 0.2 0.4 0.6 0.S 

x x 



FIGURE 5. The same as in Fig. [3|but for models (h) and (n) in Table I. Also in this case the initial state 
vectors are identical except for the values of v\. Note that in the right panel the left-propagating rarefaction 
wave covers a very small region of the flow and is closely followed by the contact discontinuity. 



of Fig. H shows that the solution to the Riemann problem still falls in the SR branch (cf. dashed 
line), with the corresponding solution at a time t > being presented in the left panel of Fig. |] 
(Note that the wave-patterns in Fig. [}]and |]both consist of a shock and a rarefaction wave, but 
have alternating initial normal velocities.). 

When nonzero tangential velocities are now considered in state 1, the right panel of Fig.^] 
shows that (^12)0 can fall in the 2R branch, hence producing a wave-pattern consisting of two 
rarefaction waves moving in opposite directions. The solution to this Riemann problem is shown 
in right panel of Fig. [| where we have chosen initial tangential velocities v\ — 0.999 and v\ = 
0.0. In this case too, it should be noted that, except for the tangential velocities, the solutions in 
the two panels of Fig. | have the same initial state-vectors but different intermediate ones. 

In analogy with Fig. ij we have collected in Fig. ||a sequence of solutions for the pressure in 
which v\ is gradually increased from to 0.999. Here too, when v\ = 0, the SR wave-pattern 
is well defined and the pressure at the contact discontinuity is intermediate between p\ and p 2 - 
However, as v\ is increased, the wave-pattern gradually changes, P3 decreases until it becomes 
smaller than p%, signalling the transition to a 2R wave-pattern. Note that while this happens, the 
region of the flow covered by the rarefaction wave becomes progressively smaller. 

In Table I we have summarized a few of the solutions shown in Figs. [| and |(| presenting numer- 
ical values for all of the relevant quantities in the Riemann problem when different combinations 
of the tangential velocities are used. 

To gain a better insight in these effects it can be instructive to consider how the velocities at 
which the various nonlinear waves propagate in the unperturbed media change when the tangen- 
tial velocities v\ and v\ are varied separately. This information is contained in Fig. ^ in which 
different curves show the behaviour of the head and tail of a left-propagating rarefaction wave 
(i.e. £y , £hdX °f tne head and tail of a right-propagating rarefaction wave (i.e. £t, and of 
a left or a right-propagating shock wave (i.e. V~, V+). The left panel of Fig. [7], in particular, 



14 



L. Rezzolla, O. Zanotti, & J.A. Pons 



1 . ! 



1 . O 




FIGURE 6. The same as in Fig. ^ but here the initial tangential velocity v{ is gradually increased from to 
0.999. The first and last solution of this sequence are also plotted in Fig. ^| 



shows the transition from a SR to a 25 wave-pattern with the dotted line marking the value of 
v 2 at which this occurs. Similarly, the right panel shows the transition from a SR to a 2R wave- 
pattern and the dotted line is again used to mark the value of v\ distinguishing the two regions 
of the parameter space. A number of interesting features can be noted and some of these were 
pointed out also by Pons et al (2000). Firstly, the speed of the head of a rarefaction wave prop- 
agating towards a region of constant tangential velocity is constant, or, stated it differently, ^ d 
does not depend on v\. Secondly, the velocity of the waves converges to zero if they propagate in 
regions with increasingly large tangential velocities. Thirdly, the values of v t at which the speeds 
of the head and tail of the rarefaction wave coincide, mark the transition from one wave-pattern 
to another and are indicated with vertical dotted lines in Fig. ^. 

As mentioned before, the appearance of these new relativistic effects is related to the behaviour 
of the function vf 2 = v'f 2 (^3) for different values of the initial tangential velocities and in partic- 
ular to how the three branches composing the curve change under variation of v\ 2 . As a result, 
the occurrence of these effects can be recast into the study of the d epend ence of(vL) 2S , (vf 2 ) 



and (vf 2 ) 2R ° n the tangential velocities. Using expressions (4.5), ( 4. 10 ), and ( 4.12 ), this depen- 
dence can be summarised as follows 



(«ia) M = («i2)2 S («a). («i2)s* = («ia)s*("i) . ("12)2* = («ia)»K. «a) > t 5 - 1 ) 



and can be best studied by considering the limits of {vf 2 ) 2S , {vf 2 ) SR and (vf 2 ) 2R when Wi :2 
00. In the case of a 25 wave-pattern, expression (4.5) simply indicates that 



lim (v^ 2 ) 2 

W2— >oc 



. 



(5.2) 



This result is also shown in the left panel of Fig. |^, where the right solid dot converges to zero 
as W 2 — > 00, while the left one does not vary. The limit (5.2) can also be used to deduce that for 
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Model 


v? 


vE 


v\ 


v\ 


p* 


«; 


p3 


P3> 


wave-pattern 


(a) 


0.5 


0.0 


0.0 


0.000 


0.597 


0.640 


0.734 


0.342 


SR 


(b) 


0.5 


0.0 


0.0 


0.300 


0.621 


0.631 


0.751 


0.349 


SR 


(c) 


0.5 


0.0 


0.0 


0.500 


0.673 


0.611 


0.788 


0.364 


SR 


(d) 


0.5 


0.0 


0.0 


0.700 


0.787 


0.570 


0.866 


0.394 


SR 


(e) 


0.5 


0.0 


0.0 


0.900 


1.150 


0.455 


1.088 


0.474 


2S 


(f) 


0.5 


0.0 


0.0 


0.990 


2.199 


0.212 


1.593 


0.647 


2S 


(g) 


0.5 


0.0 


0.0 


0.999 


3.011 


0.078 


1.905 


0.750 


2S 


(h) 


0.0 


0.5 


0.000 


0.0 


0.154 


0.620 


0.326 


0.162 


SR 


(i) 


0.0 


0.5 


0.300 


0.0 


0.139 


0.594 


0.306 


0.152 


SR 


0) 


0.0 


0.5 


0.500 


0.0 


0.115 


0.542 


0.274 


0.136 


SR 


(k) 


0.0 


0.5 


0.700 


0.0 


0.085 


0.450 


0.228 


0.113 


2R 


(1) 


0.0 


0.5 


0.900 


0.0 


0.051 


0.280 


0.168 


0.084 


2R 


(m) 


0.0 


0.5 


0.990 


0.0 


0.031 


0.095 


0.123 


0.061 


2R 


(n) 


0.0 


0.5 


0.999 


0.0 


0.026 


0.031 


0.110 


0.052 


2R 



TABLE 1. Solution of the modified Sod's problem at t = 0.4. All models refer to an ideal EOS with 
7 = 5/3 and share the same values of pressure and rest-mass density: pi = 1.0, pi = 1.0, p2 = 0.1, 
P2 — 0.125. The only differences present in the problems considered are in the normal relative velocity 
and in the tangential velocities. These quantities are reported in the first three columns, while the remaining 
ones show a few relevant quantities of the solution in the newly formed region as well as the wave pattern 
produced. 
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FIGURE 7. Velocities of the various nonlinear waves when the tangential velocities v\ and v\ are varied 
separately. The initial conditions are those of Sod's problem and the different curves refer to the head 
and tail of a left-propagating rarefaction wave (i.e. £y , £^" d ), to the head and tail of a right-propagating 
rarefaction wave (i.e. , £hd)' ar, d to a ^ l or a right-propagating shock wave (i.e. , V^). 



any (vf 2 )v > 0, there exists a value W$ of W2 such that 

« 2 )o > (vf 2 ) 2S for W 2 >W 2 . 



(5.3) 



A direct consequence of (5.3) is that given a Riemann problem having initial state vectors with 
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positive relative normal-velocity and producing a SR wave-pattern, it is always possible to trans- 
form it into a 25 wave-pattern by increasing the value of the initial tangential velocity in the state 
of initial lower pressure. 

In the case of a SR wave-pattern we refer to (4.10) to notice that in the limit of W% — ► oo the 
integrand vanishes (Ai — > oo) and therefore: 



Jim (vf 2 ) SR = 



(5.4) 



As for the previous one, the limit (5.4) can be deduced from the right panel of Fig. |2|, where the 
left solid dot converges to zero as W% — > oo, while the right one does not vary. Also in this case 
the limit (5.4) can be used to conclude that for any (vf 2 )o < 0, there exists a value W\ of W\ 
such that 

(«i 2 )o < fc > r W 1 >W 1 , (5.5) 

and therefore causing an initial SR wave-pattern solution to become a 2R one as a consequence 
of an increased tangential velocity in the state of initial higher pressure]]] 

Overall, expressions ( p".2| ) and (5.4) indicate that for tangential velocities assuming increas- 
ingly larger values, the SR branch of the vf 2 curve spans a progressively smaller interval of 
relative normal-velocities. When the tangential velocities reach their asymptotic values, the SR 
branch reduces to a point. In practice, therefore, the main effect introduced by relativistic tan- 
gential velocities in a Riemann problem is that of disfavouring the occurrence of a wave-pattern 
consisting of a shock and a rarefaction wave. 

For completeness we also report the limit of the relative normal-velocity marking the branch 
of two rarefaction waves separated by vacuum. In this case, the limit is taken for both W\ and 
W% tending to infinity, and using (4.13)-(4.14) yields 



T 1™ K 2 ) 2 h = 

Wi 2^00 



(5.6) 



Because the effects discussed in this Section have a purely special relativistic origin they might 
conflict with our physical intuition, especially when the latter is based on the knowledge of 
the Riemann problem in Newtonian hydrodynamics. However, the behaviour reported here is 
typical of those special relativistic phenomena involving Lorentz factors including also tangential 
velocities. A useful example in this respect is offered by the relativistic transverse-Doppler effect, 
in which the wavelength of a photon received from a source moving at relativistic speeds changes 
also if the source has a velocity component orthogonal to the direction of emission of the photon 
(Rindler, 1980). In this case too, a Lorentz factor including the transverse velocity is responsible 
for the effect. 

Finally, it should be pointed out that there exists a set of initial conditions for which these 
new relativistic effects cannot occur. These initial conditions are those in which vf = v 2 as in 
the classic "shock-tube" problem, where vf = = v 2 . In this case, in fact, (^12)0 = and, 
because of the limits (5.2) and (5.4), the solution of the Riemann problem will be given by a 
wave-pattern consisting of a shock and a rarefaction wave, independently of the values of the 
tangential velocities. 



6. Conclusions 

We have shown that and efficient solution of the exact Riemann problem in multidimensional 
relativistic flows can be obtained after exploiting the properties of the invariant expression for the 
relative normal-velocity between the two initial states. The new procedure proposed here is the 
natural extension of a similar method presented for the exact solution of the Riemann problem 

f Note that it is not possible to provide analytic expressions for W2 nor for Wi. 
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in one-dimensional relativistic flows (Rezzolla & Zanotti, 2001). Using information contained in 
the initial state vectors, this approach predicts the wave-pattern that will be produced in the Rie- 
mann problem, determines the set of equations to be solved and brackets the interval in pressure 
where the solution should be sought. Because logically straightforward, this approach results into 
an algorithm which is very easy to implement numerically and work is now in progress to assess 
the computational speedup in multidimensional codes. 

An important aspect of this strategy is that it naturally points out relativistic effects that can 
take place whenever the initial relative velocity normal to the initial surface of discontinuity is 
nonzero. When this is the case, in fact, the tangential velocities can affect the solution of the 
Riemann problem and cause a transition from one wave-pattern to another one. More specifi- 
cally, by varying the tangential velocities on either side of the initial discontinuity while keeping 
the remaining state vectors unchanged, the nonlinear waves involved in the solution Riemann 
problem can change from rarefaction waves to shock waves and vice-versa. These effects have a 
purely relativistic nature, do not have a Newtonian counterpart and could be relevant in several 
astrophysical scenarios, such as those involving relativistic jets or 7— ray bursts, in which nonlin- 
ear hydrodynamical waves with large Lorentz factors and complex multidimensional flows are 
expected (Blandford 2002, Meszaros 2002). 

As a final remark it is worth pointing out that while the content of this paper is focussed on 
special relativistic hydrodynamics and fiat spacetimes, the local Lorentz invariance allows to 
extend the results discussed here also to curved spacetimes and general relativistic numerical 
calculations (Pons et al 1998, and Font 2000 for a review). 



Appendix A: Monotonicity of the relative velocity as function of 

This Appendix is devoted to the proof that vf 2 is a monotonic function of p*; as mentioned in 
the main text, this is an important property and the basis of our approach. 

With our choice of considering the initial left state as the one with highest pressure, the proof 
of monotonicity will be obtained if we show that vf 2 is a monotonically increasing function of 
p*. Indicating then with a ' the first derivative with respect to p* and dropping the upper index x 
in the notati on f or the normal velocities, it is straightforward to derive that the first derivative of 



expression ( |3.6[ ) is given by 



i 1 - vie) ~ v' 2 Al ~ vie) 



12 — 7i Y2 

(1 - Vxfi V2,C) 



(1) 



A rapid look at expression ([]]) suggests that the proof that vf 2 is monotonically increasing will 
follo w if it can be shown that v[ c and — v' 2 c are both positive. On the other hand, using equations 
( |3.4| ) and (3^5), the proof of the monotonicity will follow from showing that (wf )' < and 



( v 3>)' > 0. While these inequalities must hold irrespectively of the nonlinear wave considered, 
the proofs will be different for the different waves considered. 

In the case when a rarefaction wave is present, the proof is indeed straightforward. According 



to ( 3.16 ) and ( 3.19 ), in fact, (v x )' across the rarefaction wave is negative when the rarefaction 
wave propagates towards the left [implying that (wf )' is negative] and it is positive when the 
rarefaction wave propagates towards the right [implying that (v^,)' is positive]. 

If a shock wave is present, on the other hand, a proof for the most general case and in terms of 
simple algebraic relations and cannot be given. On the other hand, a rather simple analytic proof 
can be found in the simpler case in which v'f = v 2 = 0; while this is certainly not the most gen- 
eral case, numerical calculations have shown that the result holds in general. Consider therefore 
a shock wave propagating towards the leflfjj after lengthy but straightforward calculations it is 

t A similar analysis can be repeated for the right-propagating shock wave 
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possible to show that 

(«!)' = 



H^Vs - vf)(l - V s vf) - ffxAKl - {vf) 2 W - K'(Ap) 2 



[Hi(V s 



A P V S ] 



(2) 



where we have set Hi = h\p\W 2 and Ap = p — pi > 0, and where p is the pressure behind the 
shock. If we now impose that vf — 0, we can write the derivative of (3.8) as 

V' s J' P \wl 



v s j j 2 + p \w 2 ■ 

Inserting ([3]) in (|2|) we can conclude that (uf )' is negative if and only if 
J{p\Wf + J 2 )h x - Pl ApJ' { Pl h x Wl + Ap) < 
Using ( 3.13 ) to calculate J', one finds that (|J) can be written as 



giCgi-Ap) 
j 2 



- 21l! < 



1 pxiHx+Ap) 



7(7- !) 



- - 7(7 - 

e 



2) 



(3) 



(4) 



(5) 



where p and e are the rest-mass density and the specific internal energy behind the shock front. 
Because the right hand side of ([5]) is always positive for any 7 < 2, the condition for monotonic- 
ity ([5]) will be satisfied if its left hand side is negative, i.e. if 

Ap - H x _ 



— J < -Pi- 



2hi 



(6) 



At this point the proof can be continued graphically and making use of the Taub adiabat. In the 
plane (h/p,p), in fact, the Taub adiabat (3.14) selects the points solutions of the hydrodynamical 
equations across a shock wave, therefore connecting the state ahead of the shock front with the 
one behind it. In Fig. |[ this curve is indicated with a solid line and we have indicated with the 
points A and B the states ahead (region 1) and behind (region 3) the shock front. Once an initial 
state A has been chosen, the mass flux will determine the point B of the Taub adiabat solution 
of the Rankine-Hugoniot relations. Because of this, the slope of the chord connecting the point 
A with B (shown as a dotted line in Fig. ^ is equal to — J 2 . Indicated with a dashed line, Fig. [| 
also shows the equivalent of the Taub adiabat passing through the state A but having mass flux 
equal to a 1 / 2 , i.e. 

_hi 
,P Pi. 

The point B' on such a curve represents the state behind the shock wave and, as it is clearly 
shown in Fig. |i| the slope of the chord AB' is always larger than the corresponding slope for 
the chord AB, thus stating that the condition (|^) is indeed verified and that (ff/)' is therefore 
positive. 



P : 



Pi 



(7) 



Appendix B: An explicit expression for V s 



In this Appendix we provide an explicit expression of the velocity of the shock wave propa- 
gating towards the right in the limit of p% — > p\ and when the fluid is ideal. This quantity, which 
is necessary to calcu late the limiting relative velocity («i2) 2S m equation (4.5), can be easily 
computed as [cf. eq. (3.8)] 
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FIGURE 8. The solid curve represents the Taub adiabat and the points A and B on it the states ahead and 
behind a left-propagating shock wave. See the main text for a discussion. 



where the mass flux J23' is given by [cf. eq. (3.13)] 



I 2 - 

J 23' — 



7 



Pi ~P2 



7-1/ hvQiv - 1)M - h 2 (h 2 - l)/p 2 



k3 ' 2[(j-l)p 2 + Pl ] 
The quantity V in the root dTol) is just a shorthand for 



V = 1 - 4 7 p! 



(7 - l)p 2 +Pl 
h-l) 2 ( Pl ~P2) 2 



h 2 (p2 - Pi) 



f>2 



-h\ 



(9) 



and where, finally, hy is the positive root of the Taub adibat (3.14) in the limit of p 3 — > pi, i.e 

(VV - 1)( 7 - l)(pi - P2 ) 



(10) 



(11) 
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